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The electronic origin of the huge magnetostructural effect in layered Fe-As compounds is eluci- 
dated using LiFeAs as a prototype. The crucial feature of these materials is the strong covalent bond- 
ing between Fe and As, which tends to suppress the exchange splitting. The bonding-antibonding 
splitting is very sensitive to the distance between Fe and As nuclei. We argue that the fragile in- 
terplay between bonding and magnetism is universal for this family of compounds. The exchange 
interaction is analyzed in real space, along with its correlation with covalency and doping. The 
range of interaction and itinerancy increase as the Fe-As distance is decreased. Superexchange 
makes a large antiferromagnetic contribution to the nearest-neighbor coupling, which develops large 
anisotropy when the local moment is not too small. This anisotropy is very sensitive to doping. 



Layered iron-pnictide compounds have recently at- 
tracted a lot of interest due to their high superconduct- 
ing transition temperaturcii Pairing is widely believed 
to be mediated by spin fluctuations in these materials;^ 
understanding of magnetic interaction is therefore of ut- 
most importance. Numerous first-principles calculations 
revealed a huge magnetostructural effect manifesting it- 
self, in particular, in the strong sensitivity of the Fe local 
moments to the Fe-As distance i?Fo-As'^ 

Tight-binding parameterizations of the band struc- 
ture indicate that Fe-As hybridization is significant,^'^ 
and it was used to explain the instability of the Fe lo- 
cal moments^i^ (a similar effect is known for zincblende 
iron pnictides^). However, ab initio results have been 
mainly described in terms of weak^ Fe-As hybridiza- 
tion comparable to oxides Strong Fe-As mixing was 
mentionedf^iii but its role in the magnetism was not ex- 
plained. Ref. 9 focused on the effect of the local moment 
on As-As bonding. In this paper we show that the cou- 
pling between the local moment and the Fe-As distance 
is controlled by strong covalent Fe-As bonding, and ana- 
lyze its effects on the exchange interaction in the (likely 
ground-state) "stripe" phase using the linear response 
technique. The salient features of chemical bonding and 
its relation to magnetism appear to be universal across 
the whole family of iron-pnictide layered materials, and 
we chose LiFeAs, which is a superconductor below 18 
Kfi^ as a representative example. 

The band structure of LiFeAs was calculated by 
Singh;— in most respects it is similar to LaFeAsO, 
BaFe2As2, and other members of this layered iron- 
pnictide family. Let us analyze the orbital content of 
the Bloch states. Fig. [1^ shows the energy bands^^ 
of non-magnetic LiFeAs calculated for the experimental 
structure. In this picture, As (4p) weight is shown in 
red color and also using line thickness; the iron (4s and 
3d) weight is shown in blue. We have verified that the 
weights of different 3d cubic harmonics on Fe sites (not 



shown) are similar to those plotted in Ref. d for LaFeAsO. 
It is seen from Fig. [T^ that As and Fe states form fully 
mixed bonding and antibonding states centered, respec- 
tively, at 3.5 eV below and at 1.5 eV above the Fermi level 
Ep. The large bonding-antibonding splitting of about 5 
eV indicates very strong covalent bonding between As p 
and Fe d states, in the sense that the splitting is large 
compared to the bare level separation. The hybridized 
bands have almost equal weights of Fe and As states, and 
the As states contribute equally to the bonding and anti- 
bonding states. This picture is somewhat different from 
that presented for LaFeAsO,— where the bands were di- 
vided in an upper group of mainly Fe 3d bands and a 
lower group of mainly pnictogen and oxygen p bands. 
The separation of As states in two subbands can only 
be explained by hybridization with Fe. Indeed, the As 
states of a fictitious system with Fe atoms removed from 
the lattice form a gapless set of bands about 4 eV wide; 
this was checked using the FLAPW method.^'* We also 
note that although the bare Fe 4s states are a few volts 
above Ep, they contribute appreciably at 5-6 eV below 
and at ~ 2.5 eV above Ep. 

Fig. [5] shows the valence charge density in the vertical 
plane cutting through the nearest Fe and As nuclei. The 
As sites are easily identified by the presence of two nodes 
in the radial Ap wavefunction. Strong covalent bonds be- 
tween Fe and As are clearly seen. The filling fraction 
of the As p states is only 39%, although the extended p 
orbital spills out of the atomic sphere (2.58 a.u.) some- 
what. A fully developed covalent bond corresponds to a 
filling fraction of 50%. Note that the Pauling electroneg- 
ativity difference between Fe and As is only 0.35, so that 
an almost non-polar bond is expected. 

For a strong covalent bond, the bonding-antibonding 
splitting is very sensitive to the overlap integrals, and 
hence to the distance between the atoms participating in 
the bond. Fig. [TJa shows the effect of moving As atoms 
closer to Fe layers so that i?Fc-As is reduced from 2.42 
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FIG. 1: Energy bands of LiFeAs in the non-magnetic state. MuUiken weights of As (4p) orbitals and Fe (3d and 4s) orbitals 
are shown by red and blue color, respectively. Bands with no contributions from these states appear black; fully hybridized 
bands appear purple. In addition, the As (4p) weight is shown by line thickness, (a) Experimental lattice parameters with 
Rfb-Ab = 2.42 A. (b) Rfb-Ab reduced to 2.33 A. (c) External potential of 5 eV applied to the As p states (see text). 




FIG. 2: Valence charge density (calculated using the FLAPW 
method) of non-magnetic LiFeAs in the plane cutting through 
Fe, As, and Li nuclei. Fe nuclei are at the corners of the plot. 
The linear color scale extends from (blue) to 0.2 a.u. (red). 

to 2.33 A. This relatively small change in i?Fo-As has a 
large effect on the band structure. An increase of the 
bonding-antibonding splitting is evident; the antibond- 
ing Fe-As states move notably upward, pulling the Fe 
partial density of states away from the Fermi level. The 
main change at the Fermi level is the upward shift of 
the dispersive band near the Z point, which contains a 
significant As Pz weight. This upward shift of the an- 
tibonding states results in the fourfold reduction of the 
local moment in the stripe^^ phase (see Table I). 

It is useful to view the Fe-As structure as being formed 
by two As semi-layers adsorbed on each side of a free- 
standing Fe monolayer (which is compressed by 6% com- 
pared to (001) layers in bulk bcc Fe). The strong hy- 
bridization of Fe with As apparent in Figs. [T^ and [2] ap- 
pears natural in view of this analogy. In fact, the exper- 
imentally measured position of the (chemically similar 
but slightly smaller) phosphorus adsorbed as a c(2 x 2) 
layer on the Fe(OOl) surface^^ translates to the Fe-P bond 
length of 2.27 A, which is strikingly close to the Fe-As 
distance in Fe-As compounds. 

Surfaces of magnetic transition metals often lose 
their magnetization under adsorption of such elements 
as O, S, H, N, P, etc. First-principles calculations 



for such surfaces often show reduced magnetic mo- 
ments in the surface layer which is strongly bound 
to the adsorbant) ^^'^^i^° this effect is also observed 
experimentallyi^i2^'2^ The magnetic "dead layer" on the 
surface appears precisely due to the chemical bonding, 
which partially removes the transition-metal 3d states 
from the Fermi level. Bonding and antibonding surface 
subbands are common in such cases.— Even if the sur- 
face retains some magnetic moment, it may be reduced 
compared to the bulk. Indeed, as the material is mag- 
netized and the Fermi level reaches the majority-spin 
antibonding states, further exchange splitting becomes 
unfavorable. All these arguments fully apply to Fe-As 
compounds. Although the surface of bcc Fe usually re- 
mains magnetic under chemisorption, the local moment 
in the "free-standing" monolayer in Fe-As compounds is 
expected to be less stable, because it is not supported by 
the magnetized bulk. 

It is reasonable to hypothesize that the emptying of 
the antibonding states is the primary factor responsible 
for the equilibrium position of As atoms above the Fe lay- 
ers. As soon as these states move above the Fermi level, 
further reduction of i?Fo-As docs not bring a large gain 
in the binding energy. The proximity of the antibond- 
ing Fe-As states to the Fermi level appears to be uni- 
versal among the iron-pnictide compounds. The analy- 
sis presented above indicates that covalent bonding com- 
petes with the tendency to form local moments in the 
Fe layer. It is therefore quite natural that the Fe mag- 
netic moment is extremely sensitive to i?Fc-As: a-s noted 
by other authors. ^'^^^ Competition between covalency and 
magnetism can also explain why band structure calcula- 
tions with optimized As positions seem to be in better 
agreement with experiment compared to those that use 
experimental As positions Indeed, while the local den- 
sity approximation or other exchange-correlation poten- 
tials used in density functional theory (DFT) may err 
in the binding energies, it is reasonable to expect that 
equilibrium As positions, whether in nature or in DFT, 
correspond to a similar balance between covalent binding 
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and the tendency to form local moments. 

We now concentrate on the magnetic properties of 
LiFeAs. The so-called stripe phase is believed to be the 
ground state for other iron-pnictide compoundsj ^^i^^ we 
therefore focus on this state. The stabilization energy at 
the experimental structural parameters is found to be 64 
meV per Fe site in FLAPW or 56 meV in LMTO-ASA; 
the local moment is 1.40 /is in FLAPW or 1.49 /is in 
LMTO-ASA. Further analysis shows that the trends in 
the exchange interaction in the stripe phase are insen- 
sitive to the particular choice of ASA parameters; the 
following results are, therefore, quite robust. 

Among the possible exchange mechanisms, superex- 
change and Fermi surface nesting between electron 
and hole sheets are usually mentioned for Fe-As 
compoundsj^!^ We calculated the pair exchange param- 
eters for LiFeAs using the linear response technique^- 
with a subsequent division by S^Sj, where Si is the total 
spin moment in the atomic sphere at site i. These param- 
eters map the energies of small deviations from the refer- 
ence state to the Heisenberg model E = Eq—'^-- JijSiSj. 
They were analyzed as a function of i?Fc-As, doping level, 
and also of the fictitious external potential V coupled to 
the occupation Up of the As p orbital (see below). 

The calculated pair exchange parameters for several 
nearest neighbors in the stripe phase are listed in Tables U 
andim The Tables also include the values Jq — J2iPoiJoi 
where pij = 1 or — 1 for parallel and antiparallel spin 
pairs, respectively; Jq is proportional to the Weiss field. 

The exchange interaction is quite long-range; some pa- 
rameters for pairs beyond those included in the table are 
comparable with, say, Jo20- The contribution of the first 
two coordination spheres to Jq declines steadily from 0.86 
at = to 0.35 at = 0.04c (TablelJ). A similar, but 
weaker trend is observed when V is increased (Table |ll]) . 
Thus, the reduction of the local moment is accompanied 
by the increase of the interaction range. Further, the ra- 
tio Jo / Jqo , where Jqo is defined similar to Jij and reflects 
the magnitude of the on-site (Hund) exchange, steadily 
grows from 0.26 at A^ = to 0.70 at A^ = 0.04c. These 
trends indicate that the degree of itinerancy significantly 
increases as i?Fc-As is decreased. 

When the local moment is small, the exchange param- 
eters are isotropic, as required by symmetry. For larger 
moments there is a large anisotropy, in particular for 
nearest neighbors (NN). For all Fe-As distances (Table |T| 
the NN and next-nearest (NNN) exchange parameters are 
antifcrromagnetic (AFM); the exchange along the stripe 
is thus frustrated. The ratio Jno/Jioo decreases from 
0.68 at A^ = (see TablelJ) to 0.39 at A^ = 0.04c; this 
trend tends to make the stripe phase less stable. 

We now discuss the dependence of exchange coupling 
on the external potential V coupled to the occupation of 
the As p states. This is done by adding an additional 
term Vup to the Hamiltonian; DFT self-consistency is 
achieved for each value of V. Adding V results in the 
raising of both bonding and antibonding Fe-As states, 
so that the antibonding states are moved away from the 



Fermi leveli^s This change is expected to have a strong ef- 
fect on the magnitude of the superexchange interaction!^ 
Fig. [TJ; shows the energy bands obtained with V = 5 eV. 
The band structure close to the Fermi level is very similar 
to V — (Fig. [T]). Therefore, direct exchange (includ- 
ing the "nesting" contribution) should be insensitive to 
V. Thus, the analysis of trends associated with the As 
p level shift will allow us to disentangle direct metallic 
exchange from superexchange contributions. 

Table |TT] shows that the increase of V leads to the de- 
crease and eventual sign change of the NN coupling to 
ferromagnetic. The NNN exchange, however, is almost 
unaffected. This result suggests that NN exchange has a 
large AFM contribution from superexchange, while Jno 
is primarily due to direct metallic exchange including the 
nesting effects. Moreover, the direct exchange contribu- 
tion to the NN exchange appears to be ferromagnetic. 

Table also shows the dependence of the exchange 
parameters on the doping level, which was changed by 
varying the charge of the Li nuclei. The local moment is 
quite insensitive to doping of ±20%, but a drastic change 
of exchange parameters is apparent. Most notably, the 
anisotropy of NN exchange interaction is reduced from 
Jino/Joio ~ 1 at 20% hole doping to 0.15 at 20% electron 
doping. Also, the ratio Jno/Jiooj which is important for 
the stability of the stripe phase, increases from about 0.5 
to 2.1. In the stripe phase the exchange splitting brings 
the Fermi level close to the antibonding states. It is likely 
that the anisotropy of NN coupling reflects the sensitiv- 
ity of superexchange to the anisotropy of the electronic 
structure. Interestingly, Jioo is much more sensitive to 
both the Fe-As distance and doping than Jqio- 

The exchange parameters between the Fe layers are 
quite small and frustrated for the stripe phase (the first 
three shown in the Tables are all AFM). These parame- 
ters likely vary between different Fe-As compounds. 

In conclusion, layered iron-pnictide compounds are 
characterized by strong covalent Fe-As bonding which is 
antagonistic to magnetism and responsible for the huge 
magnetostructural coupling. The magnetic interaction 
is long-range and shows increasing itinerancy with de- 
creasing i?Fc-As- Superexchange makes a dominant AFM 
contribution to the NN magnetic coupling, while the 
NNN coupling is mainly due to direct exchange. The 
anisotropy of NN coupling is sensitive to doping. 
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TABLE I: Pair exchange parameters in mRy as a function of i?Fo-As- is the shift of the As layers toward Fe layers given as a 
fraction of the lattice parameter c. Experimental structure corresponds to A^ =0. M is the local moment in hb- The exchange 
parameters Jr are indexed by the crystallographic indices of the connecting vector R in the simple tetragonal Fe sublattice, in 
which the x axis is aligned parallel to the stripes. A spin pair is parallel if the second index is even, and antiparallel otherwise. 



A, 


-Rfo-As 


M 


Jioo 


Joio 


JllO 


^200 


Jo20 


Jl20 


J210 


Jooi 


JlOl 


Jon 


Jo 





2.4204 


1.28 


-0.83 


-2.61 


-0.56 


-0.24 


-0.23 


0.27 


-0.036 


-0.063 


-0.024 


-0.035 


6.7 


0.01c 


2.3972 


1.05 


-1.26 


-2.96 


-0.75 


-0.19 


-0.14 


0.30 


-0.033 


-0.053 


-0.039 


-0.049 


8.1 


0.02c 


2.3739 


0.81 


-1.73 


-3.11 


-0.91 


0.03 


0.09 


0.35 


-0.009 


-0.044 


-0.044 


-0.057 


9.8 


0.03c 


2.3517 


0.57 


-2.10 


-3.10 


-1.01 


0.22 


0.34 


0.38 


0.025 


-0.032 


-0.061 


-0.052 


11.6 


0.04c 


2.3294 


0.31 


-2.43 


-2.90 


-0.94 


0.63 


0.60 


0.36 


0.18 


-0.032 


-0.016 


-0.028 


13.6 



TABLE II: Pair exchange parameters as a function of the external potential V coupled to As p occupation (see text). The last 
four lines are for V = with the doping level per f.u. listed in the first column. The notation is the same as in Table [T] 



V 


M 


Jioo 


Joio 


JllO 


J200 


J020 


J120 


J210 


Jooi 


JlOl 


Jon 


Jo 


-1 eV 


1.52 


-0.64 


-2.68 


-0.46 


-0.20 


-0.13 


0.33 


0.017 


-0.087 


-0.030 


-0.044 


6.3 





1.28 


-0.83 


-2.61 


-0.56 


-0.24 


-0.23 


0.27 


-0.036 


-0.063 


-0.024 


-0.035 


6.7 


1 eV 


1.10 


-0.70 


-2.30 


-0.65 


-0.14 


-0.22 


0.24 


-0.068 


-0.060 


-0.037 


-0.052 


7.5 


2 eV 


0.96 


-0.43 


-1.86 


-0.71 


-0.019 


-0.22 


0.21 


-0.089 


-0.058 


-0.048 


-0.073 


8.0 


3 eV 


0.84 


-0.06 


-1.33 


-0.78 


0.051 


-0.25 


0.16 


-0.13 


-0.040 


-0.060 


-0.097 


8.1 


4 eV 


0.70 


0.30 


-0.65 


-0.87 


-0.059 


-0.32 


0.033 


-0.25 


-0.0082 


-0.068 


-0.122 


7.7 


5 eV 


0.56 


0.84 


0.30 


-0.79 


-0.19 


-0.27 


-0.087 


-0.33 


-0.0013 


-0.121 


-0.154 


6.6 


-0.2e 


1.20 


-1.24 


-1.26 


-0.63 


-0.12 


-0.49 


0.12 


-0.30 


0.052 


-0.028 


-0.005 


3.0 


-O.le 


1.22 


-1.19 


-2.13 


-0.54 


-0.09 


-0.29 


0.34 


-0.19 


-0.042 


-0.034 


-0.052 


5.3 


+0.1e 


1.34 


-0.56 


-2.80 


-0.71 


-0.35 


-0.38 


0.15 


-0.02 


-0.011 


-0.0042 


-0.013 


7.0 


+0.2e 


1.40 


-0.41 


-2.65 


-0.88 


-0.41 


-0.49 


-0.038 


-0.01 


0.003 


-0.022 


-0.027 


6.7 
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